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<^> Abstract 

O ' 

An overview is given on the present status of the understanding of the origin of galactic cosmic rays. Recent measurements of 
charged cosmic rays and photons are reviewed. Their impact on the contemporary knowledge about the sources and acceleration 
mechanisms of cosmic rays and their propagation through the Galaxy is discussed. Possible reasons for the knee in the energy 
spectrum and scenarios for the end of the galactic cosmic-ray component are described. 



3 



in 

> 
— 

CO l. 
CN 
O 
l> 
O 



Key words: cosmic rays, origin, acceleration, propagation, knee, air shower 
PACS: 98.70.Sa, 96.50.sd 



Introduction 



The origin of high-energy cosmic rays is one of the open 
questions in astroparticle physics. The fully ionized atomic 
**^ nuclei reach the Earth from outside the solar system with 
O energies from the GeV range up to 10 20 eV. Most of them 
JL are assumed to originate in the Milky Way. At the high- 
5_i est energies, exceeding 10 17 eV, the particles are usually 
considered of cxtragalactic origin. This review focuses on 
C$ galactic cosmic rays. To distinguish between different mod- 
^> els of the cosmic-ray origin requires detailed measurements 
? *~j of the energy spectrum, mass composition, and arrival di- 
r> rection of charged cosmic rays. Additional and complemen- 
tary information is obtained through the measurements of 
high-energy photons up to TeV energies. 

This review starts with a short overview on detection 
methods typically applied to measure the composition of 
cosmic rays on satellites, balloons, and at ground level 
(Sect. [2]). The propagation of cosmic rays through the 
Galaxy is in the focus of Sect. [3] Recent progress con- 
cerning the sources and acceleration of the high-energy 
particles is described in Sect.[U One of the headstones in 
understanding the origin of galactic cosmic rays is to know 
the reasons for the knee in the energy spectrum. In Sect. [5] 
recent measurements of the all-particle energy spectrum, 
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the mean mass, and spectra for individual elements are 
compiled and their impact on the contemporary knowledge 
of the origin of galactic cosmic rays is discussed. The end of 
the galactic cosmic-ray spectrum and the transition to an 
extragalactic component is briefly illuminated in Sect. [5] 

2. Experimental Techniques 

The energies of cosmic rays extend from the GeV domain 
up to 10 20 eV. Within this range the particle flux decreases 
by about 30 orders of magnitude. This has implications on 
the accuracy to determine the mass composition of cosmic 
rays. At low energies the flux is large enough to build so- 
phisticated detectors with an active area of a few 100 cm 2 
to measure the abundance of individual isotopes. On the 
other hand, at the highest energies, where only a few parti- 
cles are expected per km 2 and century, huge ground based 
installations are necessary to measure secondary products 
generated by cosmic rays in the atmosphere and the (aver- 
age) mass can be estimated coarsely only. The situation is 
sketched figuratively in Fig.Q] 

At energies in the MeV range sophisticated silicon de- 
tectors, operated in outer space, like the Ulysses HET 
Simpson et all Il992h or the ACE/CRIS (0.025 m 2 sr) 



Stone et all 1998h experiments can identify individual iso- 



topes, fully characterized by simultaneous measurements 
of their energy, charge, and mass (E, Z, A). Since the par- 
ticles have to be absorbed completely in a silicon detector 
this technique works up to energies of a few GeV only. 



Preprint submitted to Elsevier 



5 February 2008 



o 



o 
w 

CD 
\— 

c 
o 

'w 
o 

Q- 

E 
o 
O 



elemental groups In A 

"V air Cerenkov air shower arrays 



elements |Z 

calorimeters, TRD: 
emulsions 



elements +Z 

B spectrometers 

isotopes (Z,A) * 

Si detectors B 




MeV 


GeV 


TeV 




EeV 


1 1 


i 


i i 


1 1 


1 1 



" 

I 

> 

CD 



12 b 

CNJ 

I 

E 

18 " 
LU 
T5 

Z 

24 ~° 

x 



10-30 



Energy 

Fig. 1. Illustrative sketch of the composition resolution achieved by 
different cosmic-ray detection techniques as function of energy. Over 
the energy range shown the flux of cosmic rays decreases by about 
30 orders of magnitude as indicated on the right-hand scale. 

In the GeV domain particles are registered with magnetic 
spectrome ters on stratospher ic balloons, like the BESS in- 
strument ( Aiima et al. , 2000h . They are the only detectors 
discussed here which are able to identify the sign of the 
particles charge. All other detectors rely on the fact that 
the specific ionization loss is dE/dx oc Z 2 , thus \Z\ is de- 
rived from the measurements. Magnet spectrometers are 
the only detectors suitable to distinguish between matter 
and antimatter as e.g. e + — e~ , p — p, or He — He. The par- 
ticle momentum is derived from the curvature of the tra- 
jectory in a magnetic field, which limits the usage of these 
detectors to energies approaching the TeV scale. 

At higher energies particles are measured with balloon 
borne i nstrum ents on circumpolar long duration flights 
dJoned . l2005h . individual elements are identified, char- 
acterized by their charge and energy. \Z\ is determined 
through dE/dx measurements. Experimentally most chal- 
lenging is the energy measurement. In calorimeters the 
particles need to be (at least partly) absorbed. The weight 
of a detector with a thickness of one hadronic interac- 
tion length (Aj) and an area of 1 m 2 amounts to about 
1 t. Due to weight limitations ac tual detectors li ke ATIC 
( Guzik et alll2004 ) or CREAM (|Seo et all liool have to 
find an optimum between detector aperture and energy 
resolution, resulting in relatively thin detectors with a 
thickness of 1.7 K (ATIC) or 0.7 \ (CREAM) only. The 
measurement of transition radiation from cosmic-ray parti- 
cles allows to build large detectors with reasonable weight. 
The largest cosmic-ray det ector ever flown on a balloon, 
the TRACER experiment ( Gahbauer et al. . 2004 ) has an 
aperture of 5 m 2 sr. During a lon g duration balloon flight 
over Antarctica (Horandel, 2006a) and another flight from 
Sweden to Alaska (jBovld . 120061 ). up to now an exposure 
of 70 m 2 srd has been reached with this experiment. With 
such exposures the energy spectra for individual elements 
can be extended to energies exceeding 10 14 eV. 

To access higher energies very large exposures are nec- 



essary. At present reached only in ground based experi- 
ments, registering extensive air showers. In the TeV regime 
(small) air showers are observed with imaging C e renko v 



T telesco pes such as th e HEGRA dAharonian et al' . 1999), 



" HES S (lHintonl.l2004). M AGIC (|Ferend . l2005h . or VERI 



TAS ([Weekes et all 120021 ) experiments. These instruments 
image the trajectory of an air shower in the sky with large 
mirrors onto a segmented camera. 

For primary particles with energies exceeding 10 eV 
the particle cascades generated in the atmosphere are large 
enough to reach ground level, where the debris of the cas- 
cade is registered in large arrays of particle detectors. Two 
types of experiments may be distinguished: installations 
measuring the longitudinal development of showers (or the 
depth of the shower maximum) in the atmosphere and ap- 
paratus measuring the density (and energy) of secondary 
particles at ground level. 

The depth of the shower maximum is measured in 
two ways. Light-integra t ing C erenkov detec tors like the 



BRANCA (|Fowler et all l200lh or TUNKA (jGress et al 
119991 ) experiments are in principle arrays of photomultiplicr 
tubes with light collection cones looking upwards in the 
night sky, measuring the lateral distribution of Cerenkov 
light at ground level. The depth of the shower maximum 
and the shower energy is derived from these observations . 



Imagin g telescopes as in the HiRes (lAbu-Zavvad et al.. . 
l2001al) or AUGER ([Abraham et all 12004 ) experiments ob- 



serve an image of the shower on the sky through measure- 
ment of fluorescence light, emitted by nitrogen molecules, 
which had been excited by air shower particles. These ex- 
periments rely on the fact that the depth of the shower 
maximum for a primary particle with mass A relates to 
the depth of the maximum for proton induced showers as 



X = X p 

max max 



AVn.4, 



(1) 



where Xg = 36 . 7 g/cm 2 i s the r adiation length in air 
( Matthewd . [20051 : iHbrandell . 12006 id) . 

An example for an air shower ar r ay is the KASCADE- 



Grande experiment ( Antoni et al. . 2003t Navarra et al 



2004), covering an area of 0.5 km . The basic idea is to 
measure the electromagnetic component in an array of 
unshielded scintillation detectors and the muons in scintil- 
lation counters shielded by a lead and iron absorber, while 
the hadronic componen t is measured in a large calorime- 
ter ( Engler et al. . Il999h . The total number of particles at 
observation level is obtained through the measurement of 
particle densities and the integratio n of the lateral den- 



sity distribution ([Antoni et alll200lh . The direction of air 



showers is reconstructed through the measurement of the 
arrival time of the shower particles in the individual de- 
tectors. Measuring the electron-to-muon ratio in showers, 
the mass of the prima ry can be inferr ed. A Heftier model 
of hadronic showers ( Hdrandel 2006bl ) yields the relation 
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Fig. 2. Ex posure of cosm i c-ray experiments as function of energy 
for CRIS lYanasak et all |20Q]J| TRACER antarctic and Sweden 
LDB flights llHorandei ~ l2006al ; iBoylel . l200d), MACRO /E AS-TOP 
' Aglietta et all l2004bll BLANCA dFowler et all l200lll Tibet 



..Amenomori et all" 2003). E AS-TOP llAglietta et all 1 1999)) . KAS- 
CADE iAntoni et all I2005F ). KASCADE-Grande (estimated 3 yr) 
Navarra et alll2004). AGASA, and AUGER south estimated 10 yr 



Abraham et all l2004h . The grey lines are oc E 2 . 



or lg 



Ne 

N„ 



C — 0.065 In A 



This illustrates the sensitivity of air shower experiments 
to In A. To measure the composition with a resolution of 1 
unit in In A the shower maximum has to be measured to an 
accuracy of about 37 g/cm 2 (see ([T])) or the N e /N^ ratio 
has to be determined with an relative error around 16% (see 
(J5J). Due to the large intrinsic fluctuations in air showers, 
with existing experiments at most groups of elements can 
be reconstructed with A In A ss 0.8 — 1. 

The exposure achieved by several experiments is shown in 
Fig.O The integral cosmic-ray flux roughly depends on en- 
ergy as oc E~ 2 . Hence, the grey lines (oc E~ 2 ) serve to com- 
pare the exposure of detectors in different energy regimes. 
It can be recognized that the different experiments, despite 
of their huge differences in physical size, ranging from the 
~ 300 cm 2 CRIS space experiment to the w 3600 km 2 
AUGER air shower experiment, have about equal effective 
sizes when the decreasing cosmic-ray flux is taken into ac- 
count. 

3. Propagation 

The propagation of cosmic rays in the Galaxy can be de- 
scribed quantitatively, taking into account energy loss, nu- 
clear interactions, radioactive decay, spallation production, 
and production by radioactive decay, using the equation 
(|Simpsonl . [l983h 



dx 



d_ 
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dE J 

dx 1 



N A a\Ji 
A 



Ji 



Ji 



(3) 



Ji(x) is the flux of species i after propagating through an 
amount of matter x [g/cm 2 ] subject to the reduction of 
Ji(0), that represents the source term. The other parame- 
ters are: dE / dx, the rate of ionization energy loss; erf, the 
total inelastic cross section for species z; Tj, the mean life- 
time at rest for radioactive decay for species i; aij , the cross 
section for production of species i from fragmentation of 
species j; Ty, the mean lifetime at rest for decay of species 
j into species i; Na, Avogadro's Number; A, the mean 
atomic weight of the interstellar gas; and 7, the Lorentz fac- 
tor. The exact spatial distribution of the matter traversed 
becomes important only when considering radioactive iso- 
topes, since the distance traversed in a time T is p[3cT . 

Abundances of stable secondary nuclei (produced during 
propagation) indicate the average amount of interstellar 
medium traversed during propagation before escape. Be- 
cause secondary radioactive cosmic-ray species will be cre- 
ated and decay only during transport, their steady-state 
abundances are sensitive to the confinement time if their 
mean lifetimes are comparable or shorter than this time. 

At low energies around 1 GeV/n the abundance of indi- 
vidual elements in cosmic rays has been measured. A com- 
pilation of data from various experiments is presented in 
Fig. [31 The data are compared to the abundance in the so- 
lar system. Overall, a good agreement between cosmic rays 
and solar system matter can be stated. However, there are 
distinguished differences as well, giving hints to the accel- 
eration and propagation processes of cosmic rays. 

The elements lithium, beryllium, and boron, the elements 
below iron (Z = 26), and the elements below the lead group 
(Z = 82) are significantly more abundant in cosmic rays 
as compared to the solar system composition (Fig. [3]). This 
is attributed to the fact, that these elements are produced 
in spallation processes of heavier elements (fourth term in 
(|3J|) during the propagation through the Galaxy. 

Measured ratios of secondary to primary nuclei, namely 
boron/c arbon and (scandium+ vanadium)/iron, are shown 
in Fig.Hl( Yanasak et al.l . l200lh . These ratios can be used to 
estimate the matter traverse d by cosmic rays in the Galaxy 
( Garzia-Munoz et al. . 19871) . The data can be successfully 



described using a Leaky Box model, assuming the escape 
path length for particles with rigidity R and velocity (3 = 

v/c as 



26.7/3 g/cm 2 



\1.0 GV j ~ \ IA GV ) 



-1.4 



(4) 



with 8 = 0.58 (|Yanasaket all 12001]) . The lines in Fig.g] 
correspond to this function with a path length at 2 GeV/n 
around 13 g/cm 2 decreasing to w 2 g/cm 2 at 100 GeV/n. 

However, at high energies a path length according to 
(0| decreases as X esc oc E~ s and reaches extremely small 
values. They should result in large anisotropies of the ar- 
rival direction of cosmic rays, not observed by experiments 
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Fig. 4. Abundance ratio of boron to carbon and scandium + vana- 
dium to iron in cosmic rays as fun ction of energy as measured by the 
ACE /CRIS iYanasak et all l200lll and HEAC-3 llEneelmann et all 
Il990h experiments. The curves correspond to a Leaky Box model J4J. 



(jPtuskinl . 119971 : iHorandel et all 120071 ). To sustain a suit 



able path length at high energies a residual path length 
model has been proposed, assuming the relation A rp = [6.0- 
(R/10 GVr - 6 + 0.013] g/cm 2 for the escape path length 
( Swordvl . ll995( ). Recent measurements of the TRACER ex- 
periment yi eld an upper limit for the constant term of 
0.15 g/cm 2 (|Muller et all l2005h . 



The spallation processes during the cosmic-ray prop- 
agation influence also the shape of the spectra. Usually, 
it is assumed, that the energy spectra of all elements 
have the same slope, i.e. the same spectral index at the 
source. Taking the energy dependence of the spallation 
cross sections and the pathlength of the particles in the 
Galaxy into account, it is found that at Earth the spectra 
of heav y nuclei should be nat ter as compared to light el- 
ements (jHorandel et all 120071 ) . For example, the spectral 
index for iron nuclei should be about 0.13 smaller than 
the corresponding value for protons. Indeed, direct mea- 
surements indicate that the spectra of heavy elements ar e 
flatter as compared to light elements (jHorandell . 12003a ). 
e.g. the values for protons 7 p = 2.71 and iron jp e — 2.59 
differ as expected. This effect could be of importance for 
the propagation of ultra-heavy nuclei and their possible 
contribution to the all-particl e spectrum to expla in the 
second knee around 400 PeV (jHorandel et alJ . l2007ft . 

During the propagation also radioactive secondary nu- 
clei are produced (fifth term in ([3])). Their abundance can 
be used to estimate the spatial distribution of the matter 
traversed or the cosmic-ray confinement time in the Galaxy 
if the half-life time is compa r able to the confinement time 
(iGarcia-Munoz et al. . 19751 : Simpson and Garcia-Muno3 . 

Measuring the abundance of the isotopes 10 Be 
(r 1/2 = 1.5M0 6 yr), 26 Al (r 1/2 = 8.73-10 5 yr), 36 C1 (r 1/2 = 
3.07 ■ 10 5 yr), and 54 Mg (r 1/2 = (6.3 ± 1.7) • 10 5 yr) with 
the ACRE/CRIS experiment, the propagation time in the 
Galaxy has been determined to be r esc = (15.6±1.6)-10 6 yr 



4 



(jYanasaket all 120011) . 

In the Leaky Box model, the product pisMT esc is pro- 
portional to the escape path length A esc = T esc pisM/3c. 
Knowing X esc and r esc , thus allows to determine the av- 
erage density of the interstellar matter. Measurements of 
the CRIS experiment yield an average hydrogen number 
density in the interstellar matter Pism = 0.34 ± 0.04 H 
atoms/cm 3 . A comparison to the average density in the 
galactic disc pgd ~ 1/cm 3 indicates that cosmic rays spend 
most of their propagation time outside the galactic disc, 
diffusing in the galactic halo. The heigh t of the diffusion re- 
gion in to the halo has been estimated bv lMolnar and Simon 
( 2003f) with measurements of the 10 Be/ 9 Be-ratio by the 
ISOMAX experiment to be a few (« 1 — 4 kpc). Direct 
evidence of cosmic rays propagating in the galactic halo is 
obtained from the observation of the diffuse 7-ray back- 
ground, ext ending well above the dis c, by the EGRET 
experiment ( Strong and Mattox . 19961 ). The 7-ray energy 
spectrum exhibits a structure in the GeV region, which is 
interpreted as indication for the interaction of propagating 
cosmic rays with interstellar ma tter (CR+ISM— > 7r° — > 77) 
( Strong and Moskalenkolll999l ). 

At energies in the air shower domain measurements of 
the boron-to-carbon ratio (or other ratios of secondary to 
primary elements) are not feasible any more. However, in 
this energy region propagation models can be constrained 
by information about the arrival direction of cosmic rays. 

An anisotropy is expected due to the motion of the ob- 
server (at the Earth) relative to th e cosmic-ray gas, known 
as th e Compton Getting effect (|Compton and Getting . 
1935). Such an anisotropy, caused by the orbital motion 
of the Earth around the Sun has been observed for cosmic 



rays with energies of a fe w TeV (jCutler and Grooml . 11986; 



Amenomori et all l2004al ) 



The Super-Kamiokande experiment investigated large- 
scale anisotropics in the arr ival direction of cosmic rays 
with energies around 10 TeV ( Ovama . 20061 : Guillian et al. , 
2005). A 3er excess ("Taurus excess") is found with an 
amplitude of 1.04 ± 0.20 • 10~ 3 at right ascension a = 
75° ± 7° and declination S = -5° ± 9°. On the other hand, 
a deficit ("Virgo deficit") is found with an amplitude of 
-(0.94±0.14)-10~ 3 ata = 205°±7° and S = 5°±10°. The 
anisotropy observed is compatible with a Compton Get- 
ting effect caused by a velocity of about 50 km/s relative 
to the rest frame. This velocity is smaller than the rota- 
tion speed of the solar system around the galactic center 
(« 220 km/s). This implies that the rest frame of cosmic 
rays (presumably the galactic magnetic fields) is corotating 
with the Galaxy. 

The Tibet exper iment reports anis otropies in the same 
regions on the sky (| Amenomoril . 120061 ) . For energies below 
12 TeV the large-scale anisotropies show little depen- 
dence on energy, whereas above this energy anisotropies 
fade away, consistent with measurements of the KAS- 
CADE experiment in the energy range from 0.7 to 6 PeV 



tic center would cause an excess at a « 315°, 6 = 40° and 
a minimum at a — 135°, i5 = —49° with an amplitude of 
0.35%. However, the measurements at 300 TeV yield an 
anisotropy amplitude of 0.03% ± 0.03%, consistent with an 
isotropic cosmic-ray intensity. Hence, a galactic Compton 
Getting effect can be excluded with a confidence level of 
about 5a. This implies, similar to the result of the Super- 
Kamiokande experiment, that galactic cosmic rays corotate 
with the local galactic magnetic field environment. 

Upper limits of anisotropy amplitudes in the PeV re- 
gion are compatible with the anisotropies expected from 
diffusion models of cosmic-ray propagation in the Galaxy, 
while Leaky Bo x models pred i ct a too high level of 
anisotrop y (e.g 



Horandel et all 120071 ). In particular, a 



model by ICandia et al.l (|2003l ). based on an approach of 
Ptuskin. et al.l ( 19931) . is w ell compatible with the measure- 
ments (jMaier et al.L [2005h . 



4. Sources and Acceleration of Cosmic Rays 

The energy density of cosmic rays amounts to about 
p cr w 1 eV/cm 3 . This value is comparable to the energy 
density of the visible star light p s i w 0.3 eV/cm 3 , the 
galactic magnetic fields B 2 /2po w 0.25 eV/cm 3 , or the mi- 
crowave background p^K ~ 0.25 eV/cm 3 . The power re- 
quired to sustain a constant cosmic-ray intensity can be es- 
timated as L cr = p cr V/T esc m 10 41 erg/s, where r esc is the 
residence time of cosmic rays in a volume V (the Galaxy, or 
the galactic halo). With a rate of about three supernovae 
per century in a typical Galaxy the energy required could be 
provided by a small fraction (rj 10%) of the kinetic energy 
released in supernova explos i ons. T his has been realized al- 
ready bv lBaade and Zwickv (1934). The actual mechanism 
of acceleration remained mysterious until Fermi ( 1949) pro- 
posed a process that involved interaction of particles with 
large-scale magnetic fields in the Galaxy. Eventually, this 
lead to the currently accepted model of cosmic-ray accel- 
eration by a first-order Fermi mechanism that operates in 
strong shock fronts which are powered by the explosions 
and propagate from t he supernova remnant (SNR) into the 



interstellar medium (lAxford et al 



(jAntoni et all l2004af) . A Compton Getting effect caused 



by the orbital motion of the solar system around the galac- 



.19771 iKrvmskvl 11977 
19781 ). 



Belll978tlBfanford and Ostriker! 

Diffusive, first-order shock acceleration works by virtue 
of the fact that particles gain an amount of energy AE oc 
E at each cycle, when a cycle consists of a particle pass- 
ing from the upstream (unshocked) region to the down- 
stream region and back. At each cycle, there is a probabil- 
ity that the particle is lost downstream and does not return 
to the shock. Higher energy particles are those that remain 
longer in the vicinity of the shock and have enough time 
to achieve the high energy. After a time T the maximum 
energy attained is E max ~ Zef3 s BTV s , where (3 S = V s /c 
is the velocity of the shock. This results in an upper limit, 
assuming a minimal diffusion length equal to the Larmor 
radius of a particle of charge Ze in the magnetic fields B 
behind and ahead of the shock. Using typical values of type 
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measurements of the abu ndances of 59 Ni and 59 Co by the 
ACE/CRIS experiment (jWiedenbeck et all [l999h . These 



2 8 

Solar System Abundance ( Si=l) 

Fig. 5. Comparison of derived galactic cosmic-ray (GCR) source 
abundances of refractory nuclides with solar-system abundances 
according to measureme nts with ACE/CRIS normalized to 28 Si 
ijWiedenbeck et allEooS l. 



II supernovae exploding in a n average interstellar medium 
yields E max w Z ■ 100 TeV (jLagage and Cesarskvl Il983l ). 
More recent estimates give a maximum energy up to one 
order of magnitude l arger for s ome t ypes of supernovae 
E mnf r , « Z ■ 5 PeV ( Berezhko , 19961 Kobavakawa et all 
2002tlSveshnikovaL f2003) . As the maximum energy depends 
on the charge Z 1 heavier nuclei (with larger Z) can be accel- 
erated to higher energies. This leads to consecutive cut-offs 
of the energy spectra for individual elements proportional 
to their charge Z, starting with the proton component. 

The overall similarity between cosmic rays and matter 
in the solar system, as already seen in Fig.[3j indicates that 
cosmic rays are "regular matter", but accelerated to ex- 
treme energies. 

Detailed information on the composition at the source 
can be obtained from measurements of the abundance of 
refractory nuclei. They appear to have undergone minimal 
elemental fractionation relative to one another during the 
propagation process. The derived abundance at the source 
is pr esented in Fig.[5lversus th e abundance in the solar sys- 



tem (jWiedenbeck et all 120031 ) . The two samples exhibit an 



extreme similarity over a wide range. Of the 18 nuclides 
included in this comparison, only 58 Fe is found to have an 
abundance relative to 28 Si that differs by more than a fac- 
tor of 1.5 from the solar-system value. When uncertainties 
are taken into account, all of the other abundances are con- 
sistent with being within 20% of the solar values. This in- 
dicates that cosmic rays are accelerated out of a sample of 
well mixed interstellar matter. 

The time between synthesis of cosmic-ray material and 
its acceleration to high energies has been estimated by 



nucleides form a parent-daughter pair in a radioactive decay 
which can occur by electron capture only. This decay can- 
not occur once the nuclei are accelerated to high energies 
and striped of their electrons. The measured abundances 
indicate that the decay of 59 Ni to 59 Co has occurred, lead- 
ing to the conclusion that a time longer than the 7.6 • 10 4 yr 
half-life time of 59 Ni elapsed before the particles were ac- 
celerated. Such long delays indicate the acceleration of old, 
stellar or interstellar material rather than fresh supernova 
ejecta. For cosmic-ray source material to have the compo- 
sition of supernova ejecta would require that these ejecta 
not undergo significant mixing with normal interstellar gas 
before ps 10 5 yr have elapsed. 

Ratios of individual cosmic-ray elemental abundances 
to the corresponding solar system abundances seem to be 
ordered by the first ionization potential (FIP) of the ele- 
ments. In the GeV energy range elements with a FIP be- 
low « 10 eV are about ten times more abundant in the 
cosmic-ray sources than are elements with a higher FIP 



(jMever et al.l . ll997l ). It has also been found that the abun- 



dance ratios of galactic cosmic rays to the solar system val- 
ues scale with the condensation temperature T c . Refractory 
elements (T c > 1250 K) are more abundant than volatile 
elements (T c < 875 K). A model which could perhaps ac- 
count for the FIP effect would have a source ejecting its 
outer envelope for a long period with a FIP selection ef- 
fect (as the sun does) followed by a supernova shock which 
sweeps up and accelerates this material. It has been pro- 
posed that the FIP effect is due to non-volatiles being ac- 
celerated as grain s, while the yolatiles are accelerated a s 
individual nuclei (jMever et all Il998t lEllison et all . Il998f ). 



Recent measurements of the TRACER experiment allow 
to investigate these effect a t higher energies, n amely at 
100 GeV/n and 1000 GeV/n ijMiiller et all . l2005r ) . In Fig.[6] 
the abundance at the cosmic-ray sources relative to the so- 
lar system abundances are shown as function of the first 
ionization potential and as function of condensation tem- 
perature. The decreasing ratio as function of FIP and the 
increase as function of condensation temperature, known 
from lower energies, is also pronounced at energies as high 
as 1000 GeV/n, as can be inferred from the figure. 

The theory of acceleration of (hadronic) cosmic rays at 
shock fronts in supernova remnants, mentioned above, is 
strongly sup ported by recent measureme nts of the HESS 
experiment (jAharonian et al. . 20041 2006h . observing TeV 
7-rays from the shell type supernova remnant RX J1713.7- 
3946, originating from a core collapse supernova of type 
Il/Ib. For the first time, a SNR could be spatially resolved in 
7-rays and spectra have been derived directly at a potential 
cosmic-ray source. The measurements yield a spectral index 
7 = -2.19±0.09±0.15 for the observed 7-ray flux. The pho- 
ton energy spectrum of the supernova remnant RX J1713 is 
depicted in Fig. [7] Measurements in various energy ranges 
(ATCA at radio wavelengths, ASCA x-ray, EGRET GeV 
7-ray, CANGAROO and HESS TeV 7-ray) are compared 
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Fig. 6. Abundances of elements at the cosmic-ray sources relative to the solar system abundances as fu nction of the first io nization potential 
(left) and the condensation temperature (right) as measured by the CRN and TRACER experiments llMuller et all [2005h . 

with material in the vicinity of the source (hadr + ISM 
— > 7T° — > 77). This process is clearly dominant over elec- 
tromagnetic emission generated by inverse Compton effect 
and non-thermal Bremsstrahlung, as can be inferred from 
the figure. The results are compatible with a nonlinear ki- 
netic theory of cosmic-ray acceleration in supernova rem- 
nants and imply that this supernova remnant is an effective 
source of nuclear cosmic rays, where about 10% of the me- 
chanical explosion energy are conve rted into nuclear cos- 
mic rays (jVolk and Berezhko" , 2006). Further quantitative 
evidence for the acceleration of hadrons in supernova rem- 
nants is provided by me asurements of the HEGRA experi- 
ment dAharonianll200ll) of TeV 7-rays from the SNR Cas- 




5 10 
Photon energy Ig Ey [eV] 



Fig. 7. Spatially integrated spectral energy distribution of the super- 
nova remnant RX J1713. The solid line above 10 7 eV corresponds to 
7-ray emission from 7r°-decay, whereas the dashed and dash-dotted 
curves indicate the inver se Compton and nontherm al Bremsstrahlung 
emissions, respectively (Volk and Berezhko, 2006). 



siopeia A (jBerezhko et all I2003I ) and by measurements of 



the HESS experim ent from th e SNR "Vela Junior" (SNR 
RX J0852.0-4622) (IVdlkl . l200fih . 



to predictions of a model bv lVolk and Berez hkol (|2006l) . 

The solid line below 10 6 eV indicates synchrotron emis- 
sion from electrons ranging from radio frequencies to the 
x-ray regime. The observed synchrotron flux is used to ad- 
just parameters in the model, which in turn, is used to pre- 
dict the flux of TeV 7-rays. An important feature of the 
model is that efficient production of nuclear cosmic rays 
leads to strong modifications of the shock with large down- 
stream magnetic fields (B w 100 /J.G). Due to this field 
amplification the electrons are accelerated to lower max- 
imum energy and for the same radio/x-ray flux less elec- 
trons are needed. Consequently, the inverse Compton and 
Bremsstrahlung fluxes are relatively low only. The solid 
line above 10 6 eV reflects the spectra of decaying neu- 
tral pions, generated in interactions of accelerated hadrons 



5. The Knee in the Energy Spectrum 

5.1. All-Particle Spectrum 

As discussed above, in the standard picture of galactic 
cosmic rays, it is assumed that the particles gain energy in 
supernova re mnants and pr opagate diffusively through the 
Galaxy (e.g. iGaisseiiExM . The nuclei are accelerated up 
to a maximum energy, being proportional to their charge. 
During propagation some particles leak out of the Galaxy 
and the leakage probability is a function of the nuclear 
charge as well. This implies that the energy spectra for indi- 
vidual elements should exhibit a break (or knee structure) 
at energies being proportional to the elemental charge. 

The situation is sketched in Fig.[5] The lines indicate 
energy spectra for groups of elements following power laws 
with a break proportional to Z . In the particular example 
the cut-off energy is assumed as E% = Z-4.5 PeV, according 



7 



5 

CD 

C/3 

v 10 

CO 



LU 

o 
LU 

e 

"D 

X 
13 



10' 



+ 
_ + 
+ 

■:■ 

o 

^> 

+ 

X 



i 



AGASA 
Akeno 20 km 
Akeno 1 km 
AUGER 
BLANCA 
CASA-MIA 
DICE 

BASJE-MAS 
EAS-Top 
Fly's Eve 
GRAPES QGS 
GRAPES SIB 
Haverah Park 
Haverah Park Fe 
Haverah Park p 
HEGRA 



HiRes-l 
HiRes-ll 
Hi Res/Ml A 
KASCADE 
KASCADE 
KASCADE 
KASCADE 
KASCADE-i 



e/m QGS) 
e/m SIB) 
h/m) knee 

rande 



MSU 

■S 1 Mt. Norikura 

# SUGAR 
® Tibet ASy 



direct: 



Tibet AS7-I 
Tunka-25 



JACEE 
RUNJOB 

SOKOL 
Grigorov 




Energy E [GeV] 

Fig. 8. All-part i cle en ergy spect rum of cosmic rays, th e flux is mu ltiplied by E 3 . Resu lts from direct measurements by Grig orov et al. 
dGrigorov et all Il999l) , JACEE llAsakimori et all Il995l) , RUN JOB JPerbina et all |2005|). and SOK OL lllvanenko et all Il993» as well a s 



from the air shower experiment s AGASA llTakeda et all [2003T), Akeno 1 km 2 llNagano et aTlll984al) , and 20 k m 2 jNagano et al 
AUGE R JSommers et all 1 20051), BASJE-MAS llOgio et all 12004). B LANCA jFowler et all l200lh, CA SA-MIA jGlasmacher et al 



-A 

DICE llSwordv and Kiedal, 12000). EAS-TO P 
hadronic interactio n models I Ha^ashjetjiLk 200 




1984b), 



Eye 



1999b), 



v , Ha verah Park I Lawrence et al.. . _ 
l200Cf) . HiRes-MIA l|Abu-Zavvad et all 12001 j ). HiRe s-I llAbbasi et all 120041) . HiRes- II llAbbasi et al 



[ Corbato et all 119941), GRAP ES-3 inter preted with two 
Il99ll) and <Ave et all l2003t) . HEGRA l lArqueros et all 



KASCADE electrons and muons 
and a neural net work ana lysis 



interpreted with two hadronic interacti on models ([Antoni et all |2005|) . hadrons llHorandel et al T 

comb ining different sh ower component s j Antoni et all |2002|). KASCADE-Grande (pr eliminary) fHaungs et all l2006bf ) MSU llFomin et al.. . 
Il99ll). Mt. Norikura dlto et all ll997T). SUGAR llAnchordoqui and Goldbergt 12004), Tibet AS t llAmenomori et all l2000al) and AS7-III 



l|Amenomori et all [20031 ). Tunka-25 l|Chernov et all 120061) . and Yakutsk l|Glushkov et all l2003h. The li nes represent spectra for elemental 
groups (with nuclear charge numbers Z as indicated) according to the poly-gonato model (Horandel, 2003a). The sum of all elements (galactic) 
and a presumably extragalactic component are shown as well. The dashed line indicates the average all-particle flux at high energies. 



to the poly-gonato model ( Horandel . l2003af ) and a possible 
contribution of ultra-heavy elements [Z > 28) is shown 
as well. Summing up the fluxes of all elements yields the 
all-particle flux, indicated as "galactic" in the figure. The 
all-particle flux follows a power law with a spectral index 
7 sa —2.7 up to about 4.5 PeV and then exhibits a kink, the 
knee in the energy spectrum and continues with an index 
7 sa —3.1 at higher energies. This scenario agrees well with 
the flux measured with detectors above the atmosphere 
and with air shower experiments, shown in the figure as 
well. In the measured spectrum some structures can be 
recognized, indicating small changes in the spectral index 
7. The most important are the knee at Ek ~ 4.5 PeV, 
the second knee at E 2n d ~ 400 PeV« 92 x Ek, where the 
spectrum exhibits a second steepening to 7 rj —3.3, and 
the ankle at about 4 EeV, above this energy the spectrum 
seems to flatten again to 7 w —2.7. To understand the 
origin of these structures is expected to be a key element 
in understanding the origin of cosmic rays. 

Various scenarios to explain the knee are proposed in th e 
literature, for an overview see e.g. (|Horandeil200ll2005| ). 
The most popular approaches (maximum energy attained 
and leakage), just described above, are modeled with vary- 



ing details, resulting in slightly different spectra. But, also 
other ideas are discussed, like the rcaccclcration of particles 
in the galactic wind, the interaction of cosmic-ray parti- 
cles with dense photon fields in the vicinity of the sources, 
interactions with the neutrino background assuming mas- 
sive neutrinos, the acceleration of particles in 7-ray bursts, 
or hypothetical particle physics interactions in the atmo- 
sphere, transporting energy in unobserved channels, thus 
faking the knee feature. All scenarios result in spectra for 
individual elements with a break at energies being cither 
proportional to the nuclear charge Z or to the nuclear mass 
A which yield certain structures in the sum spectrum. To 
distinguish between the different models, measurements of 
the (average) mass of cosmic rays as function of energy are 
required, or - even better - the measurement of spectra of 
individual elements or at least elemental groups. 



5.2. Mean Mass 

Frequently, the ratio of the number of electrons and 
muons is used to determine the mass composition, see ([2]). 
Muons with an energy of several 100 MeV to 1 GeV are 
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Fig. 9. Mean logarithmic mass of cosmic rays deriv ed from the measurements of elec trons, muons, and hadrons at ground le vel. Results are 
shown from CASA-MIA llGlasmacher et all [l999aT), Chacaltaya [|Aguirre et al.1 l200oh . EAS-TOP electrons and GeV muons llAglietta et all 
l2004all. EAS-TOP/M ACRO (TeV mu onsl dAglietta et all l2004bh . GRAPES-3 data interpreted with two hadronic interaction models 
l|Havashi et all |2005|) , HEG RA CRT l|Bernlonret**an. Il998h. KA SCADE electrons and muons interpreted with two hadronic interaction 
models lAntonje^*aU, |2005|) . hadrons and muons jHorandel et al.l.ll998(). as well as an analysis combining different obscrvables with a neural 



network jAntoni et all I2002T). and SPASE/AMANDA l|Rawlins et al.U2003h. For comparison, results from direct measurements are shown 
as well from the JACEE JjACEE collaboration] , fl999t) and RUNJOB jDerbina et all [20051 ) experiments. For orientation, In A for selected 

elements is indicated on the right-hand side. 

Models: The the grey solid and dashed lines indicate spectra according to the poly-gonato model llHoran dcl. 20033). 

T op: The lines i ndicat e spectra for models explaining the knee due to the ma ximum energy attain ed dur ing the acceleration proce ss according 
to lSveshnikoval ||2003|) (— . ■ ■ -l. lBerezhko and Kscnofonto 3 {HH) (---l. lStanev et al.1 lll993h (• • iKobavakawa et al.l 1200211 (--■). 
Bottom: Th e lines indicate spectr a for models explaining the kne e as ef fect of l eakag e from the Gala xy during the propagation p rocess 
according to lHorandel et al.1 ll2007fl (— l. lQgio and Kakimotol (120031 1 (- - -). iRouleU i2004) (■■■). as well as I Volk and Zirakashvilil (12003 1 (--■). 



used b y the experim ents CASA-MIA (iGlasmacher et al 
1999al). EAS-TOP (lAdietta et al.l . |2004aD, GRAPES-3 



(jHavashi et"all l2005h . or KASCADE. To study system- 
atic effects two hadronic interaction models are used 
to interpret the data measured with GRAPES-3 and 



KASCADE (|Antoni et al.l . l2005h . High energy muons 
detected deep below rock or antarctic ice are ut ilized 
by the EAS-TOP/ MACRO (lAglietta et all . l2004bh and 
SPASE/AMANDA (jRawlins et all 120031) experiments. 
Also, the correlation between the hadronic and muonic 
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shower components has been investigated , e.g. by the KAS- 
CADE experiment ( Horandel et al. . 19981 ). The production 
height of muons has be en reconstructed by HEGRA/CRT 



I Bernlohr et al 
20031) . 



1998h and KASCADE (|Biittner et al 



Results from various experiments measuring electrons, 
muons, and hadrons at ground level are compiled in Fig. [51 
The same experimental data are presented in the upper and 
lower graphs, where they are compared to different models. 
At low energies the values for the mean logarithmic mass 
are complemented by results from direct measurements. A 
clear increase of (In A) as function of energy can be recog- 
nized. However, individual experiments exhibit systematic 
differences of about ±1 unit. Such fluctuations in (In A) are 
expected according to the simple estimate ([2]). assuming 
that the ratio of the electromagnetic and muonic shower 
components can be measured with an accuracy of the or- 
der of 16%. This uncertainty is a realistic value for the res- 
olution of air shower arrays. Of particular interest are the 
investigations of the KASCADE and GRAPES-3 experi- 
ments: interpreting the measured data with two different 
models for the interactions in the atmosphere results in a 
systematic difference in (In A) of about 0.7 to 1. 

The experimental values in Fig.[9]follow a trend predicted 
by the poly-gonato model as indicated by the grey lines in 
the figure which implies that the increase of the average 
mass with energy is compatible with subsequent breaks in 
the energy spectra of individual elements, starting with the 
lightest species. The experimental values are compared to 
astrophysical models for the origin of the knee in the figure. 
The top panel shows models which explain the knee due 
to the maximum energy achieved during the acceleration 
process. In the lower panel predictions from propagation 
models (including reacceleration during propagation) are 
compiled. Details of the individual models will be discussed 
below. In general, all these predictions indicate an increase 
of the mean mass as function of energy very similar to the 
measured values. 

Another technique to determine the mass of cosmic rays 
is to measure the average depth of the shower maximum 
(X max ). The results of several experiments are presented in 
Fig.[l0l It depicts the measured X max values as function of 
energy. Two different techniques are used in the measure- 
ments, namely the imaging technique, using telescopes to 
obtain a direct image of the shower using Cerenkov or flu- 
orescence light and the light-integrating method, in which 
the height of the shower maximum is derived from the lat- 
eral distribution of the Cerenkov light measured at ground 
level. The DICE, Fly's Eye, and HiRes experiments use the 
imaging method, while all other experiments belong to the 
second group. The Haverah Park experiment uses the rise 
time of the shower front to estimate X max . The data in 
Fig.nolshow systematic differences of rj 30 g/cm 2 at 1 PeV 
increasing to « 65 g/cm 2 close to 10 PeV. Some of the 
experimental uncertainties may be caused by changing at- 
mospheric conditions. The imaging experiments measure a 
geometrical height which has to be converted into an atmo- 
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Fig. 11. Interpreting the measured Xmax values with the two inter- 
action models, represented in Fie. 1101 by the solid and dashed lines, 
leads to differences A (In A) plotted here as function of energy. 

spheric depth. Measuring longitudinal atmospheric profiles 
during different seasons, Keilhauer et al. (2003) found that 
the atmospheric overburden for a fixed geometrical height 
(e.g. 8 km a.s.l.) varies by at least 25 g/cm 2 . The light- 
integrating technique is rather indirect, requiring an air 
shower model to convert the observed lateral distribution 
into X max . Application of different codes c an explain parts 
of the discrepancies in X max . For example. iDickinson et al 



( 2000|) found that using different interaction models in the 
MOCCA code, viz. the original and the SIBYLL model, 
the systematic error amounts to AX max w 10 g/cm 2 . 

The observed values are compared to predictions of air 
shower simulations for prima ry protons and ir on nuclei 
using the program CORSIKA faeck et al.l.ll998[) with t he 



hadro nic interaction model QGSJET 01 (jKalmvkov et al 



19971) and a modified version with lower cross sections and 
larger values for the elastici t y of the hadronic interactions 
(model 3a from iHorandell (|2003bl )). Both variations of 
QGSJET are compatible with accelerator mea surements 
within their error boundaries, for details see iHorandel 
(|2003a) . In principle, the difference between the two cases 
illustrated in the figure represents an estimate of the 
projection of the experimental errors from collider exper- 
iments on the average depth of the shower maximum in 
air showers. At 10 9 GeV the difference between the two 
scenarios for primary protons is about half the difference 
between proton and iron induced showers. This illustrates 
the significance of the uncertainties of the collider mea- 
surements for air shower observables. The lower values 
for the inelastic proton-air cross section (model 3a) are in 
good agree ment with recent measurements from the HiRes 
experiment l|BelovLl2006t Horandel, l2006d) . 

Knowing the average depth of the shower maximum for 
protons A,^ la2 . and iron nuclei X^ e ax from simulations, the 
mean logarithmic mass can be derived in the superposi- 
tion model of air showers from the measured X™f°f us- 
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Fig. 10. Averag e depth of the show er maxi mum X max as function of primary en ergy as obtained by t he experiments BLANCA iFowler et all 
l200ll). CACTI 1 [Paling et allll997l) , DICE dSwordv and KiedaLl200Cft , Fly's Eye jCorbato et al.Lll994) . Haverah Park jWatsoij.l200Ch. HEGRA 
dAraueros et al.ll2000h, H i Res/M IA jAbu-Zawad et al.Ll2000h. Mt. Lia n Wang |cha et al.Ll200lf ). SPASE/VULCAN {Dickinson et allll999T) . 



Tunka-25 llChernov et all |2006|) , and Yakutsk jKnurenko et all l200ll) . The lines indicate simulations for proton and iron induced showers 
using the CORSIKA code with the hadronic in teraction model Q GSJET 01 (solid line) and a version with lower cross sections and slightly 
increased elasticities (dashed line, model 3a in (Horandcl, 2003b)). 



ing (In A) = - X^ ax )/{X^ ax - X^ ax ) ■ lnA Fe . 

This conversion requires to chose a particular interaction 
model. The influence of different interaction models on 
the (In A) values obtained is discussed in detail elsewhere 
(jHorandell . l2003bf ). Taking the two cases shown in Fig.fTOl 
as solid and dashed lines yields differences in (In A) as de- 
picted in Fig.QTJ As expected, they grow as function of 
energy and exceed one unit in (In A) at energies above 
10 9 GeV. As mentioned above, these differences are pro- 
jections of the errors of parameters derived at accelerators 
on air shower measurements. The average depth of proton 
showers is more increased as the depth of iron induced cas- 
cades (see Fig- fTO]) . Since in the energy region between 10 7 
and 10 8 GeV the measurements indicate a trend towards a 
heavier composition, the influence of the modifications on 
the mean mass are slightly smaller in this region, resulting 
in the dip in Fig. [Til 

Applying the modification of QGSJET 01 (model 3a) to 
the data shown in Fig. [TO] yields (In A) values plotted in 
Fig. [TO] as function of energy. According to UJ), an uncer- 
tainty of about 1 Xq in the determination of X max yields 
differences in (In A) of about 1 unit. Such a scatter of the 
experimental values is visible in Fig. 1121 

Similar to the results shown in Fig. [9] also the observa- 
tions of X max yield an increase of the average mass as func- 
tion of energy. The same experimental results are shown 
in the upper and lower panel. The data are compared to 
astrophysical models proposed to explain the knee due to 



the maximum energy attained during acceleration (upper 
panel) and due to leakage from the Galaxy (lower panel). 
The general trend of the data is seen as well in the models 
shown. This conclusions, of course, depends on the inter- 
action model used to interpret the measurements. 



5.3. Spectra for Individual Elements 

A significant step forward in understanding the origin of 
cosmic rays are measurements of energy spectra for indi- 
vidual elements or at least groups of elements. Up to about 
a PeV direct measurements have been performed with in- 



strum ents above the atmosphere (e.g. Wiebel -Soth etal 
19981 ). As examples, results for primary protons, helium 



the CNO group, and iron nuclei are compiled in Figs. 1131 
and 1 141 The same experimental data are shown in both fig- 
ures. They are compared to different models proposed to 
explain the knee in the energy spectrum as discussed in de- 
tail below. 

Recently, also indirect measurements of elemental groups 
became possible. 

A special class of events, the unaccompanied hadrons 
were in vestigated by th e EAS-TOP and KASCADE exper- 
iments (|Aglietta et all 120031 : lAntoni et all l2004bl ). Simu- 
lations reveal that these events, where only one hadron is 
registered in a large calorimeter, are sensitive to the flux 
of primary protons. The derived proton fluxes agree with 
the results of direct measurements as can be inferred from 
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to lSveshnikoval ||2003|) (— . ■ • -l. lBerezhko and Kscnofonto 3 l|l999h (--A. IStanev et al.l lll993h (• ■ ■ ) . iKobavakawa et alj ]2002h (--■)■ 
Bottom: Th e lines indicate spectr a for models explaining the kne e as ef fect of l eakag e from the Gala xy during the propagation p rocess 
according to lHorandel et al.1 | |2007h ( — V lOgio and Kakimotol (120031 1 (- - -). iRouletl i2004) (• • •)> as well as lVolk and Zirakasnvtlil (12003 1 (--■). 



Figs.[T51and[Hl indicating a reasonably good understand- 
ing of the hadronic interactions in the atmosphere for en- 
ergies below 1 PeV. 

At higher energies a breakthrough has been achieved by 
the KASCADE experiment. Measuring simultaneously the 
electromagnetic and muonic component of air showers and 
unfolding the two dimensional shower size distributions, the 
energy spectra of five elemental groups have been derived 



( Ant oni et al. , 2005). n order to estimate the influence of 
the hadronic interaction models used in the simulations , 



two m odels, namely QGSJET 01 and SIBYLL (|Engel et al 



1999), have been applied to interpret the measurements. It 



turns out that the all-particle spectra obtained agree satis- 
factory well within the statistical errors. For both interpre- 
tations the flux of light elements exhibits individual knees. 
The absolute flux values differ by about a factor of two or 
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Fig. 13. Cosmic-ray energy spectra, caption 2 pages below. 
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Fig. 1131 and ll4l Cosmic-ray energy spectra for four groups of elements, from top to bottom: protons, helium, CNO group, 
and iron group. 

Prot ons: Re sults from direct me asurements a bove the atmosphe re by AM S dAlcaraz et all 2000), ATIC (Wefel et al 



IMA X dMenn et al 
JOB (Derbina et al 



l2000h . JACEE dAsakimori et all Il998ft. MASS (iBellotti et all Il999l). IPapini et all dl993h. RUN 
l2005li. R JCH-II (|Diehl et all 12003V . iRvan et all l|l972ft . lSmith et all (|l973l ). SOKOL dlvanenko et al 



2005ft . BESS dSanukiet al.L l2000h. C APRICE (iBoezio et all 12003ft. H EAT dVernois et al.l.|200lft.llchimura et alJ d 1993ft. 



" I I ' I j I I ^x^x^xxx ixx. |, | XJVJW.jy , | xvf xxxx K ixx . | ,| x. u i ixx. | 1| * I , '''''''''' 1| x » miwuw u ixx. |, 

Il993ft . IZatsepin et all ( 19931 ). an d fluxes obtained fr om indirect measur ements by KASCAP E electrons and muons for 



two had ronic interaction mod els ( Antoni et al. . 2005|) and single hadron s ( Antoni et al. . 2004bl ). EAS-TOP (electrons and 
muons) ( Navarra et all 120031) and s ingle hadr ons dAglietta et all 120 03). GRAPES-3 int erpreted with two hadronic inter- 
action mo dels dHavashi et al l l2005ft . HEGRA dAharonian et allll999TTMt. Chacaltay a dlnoue et allll997l). Mts. Fuji and 
Kanba la dHuang et all 120031 ). Tibet burst detector (HD) ( Amenomori et all lioOObft and AS7 (HP) dAmenomori et al 
2004bft . 



JAC EE dAsakimori et all Il998ft. MASS dBellotti et all 



1999), P apini et al 



JOB dDerbina et alll2005ft . lSmith et al.l (|l973h . SOKOL dlvanenko et allll993ft . lWebber et all (|l987| ) and fluxes obtaine d 



Helium: Results from direct measurements above the atmos phere by ATIC (Wefel et all [2005ft . BESS dSanuki et"al.l . 
2000ft . C APRICE dBoezio et all l2003ft. HEAT dVernois et alllioOlh. Ilchimura et alJ dl993l). I MAX dMenn et all l2000l). 



1993ft. RICH-II dDiehl et alll2003l) . RUN- 



from indirect measurements by KASCAPE electrons and muo ns for two hadronic in teraction models ( Antoni et all 2005ft 



GRAPES-3 interpreted with two hadro nic interaction models dHavashi et all 12005). Mts. Fuji and Kanbala dHuang et al 



2003), and Tibet burst detector (HP) dAmenomori et alll2000bft 



CNO group : Resu lts from direct measurements abov e the atmosphere by ATIC (C+O) dCherrv ) . | 2006ft . CRN (C+O) 
dMiiller et allll99lft TRACER (O) dMiiller et alll2005j ) JAC EE (CNO) d JACEE collaborationl . ll999l ). RUN JOB (CNO) 



dPerbina et al.l . l2005ft . SOKO L (CNO) dlvanenko et al 



1993ft . and fluxes obta i ned fr om indirect measurements by KAS- 



CAPE electrons and muons ( Antoni et all 20051) . GRAPES-3 dHavashi et al. . 2.005), the latter two give results for two 



hadronic interaction models, and EAS-TOP ( Navarra et al. . 2003ft . 
Iron: Re sults from direct measurements above t h e atmosphere by ATIC ( Cherry . 2006ft. CRN dMiiller et ah . 1991 ). 
HEAO-3 dEngelmann et allll985l).|juliussonl dl974l).lMinagawa! dl98lft, TRACER dMiiller et al l|2005ft (sin gle element res- 



olutio n) and lHarevama et alJ d 1999ft . I lchimura et al.l (|l993h . JACEE ( Asakimori et al J,|l995|), RUN JOB dPerbina et al" 
2005), SOKOL dlvanenko et al.L Il993h (iron group), as well as fluxes from indirect measurements (iron group) by EAS - 
TOP dNavarraet_a lll2003l), KA SCAPE electrons and muons d Antoni et alll2005ft . GRAPES-3 dHavashi et alll2005ft . 
and HESS direct Cerenkov light ( Aharonian et all 12007 ). The latter three experiments give results according to interpre- 
tations with two hadronic interaction models. 

Models: The the grey solid lines indicate spectra according to the poly-gonato model dHor andell . l2003ah . 
Fig. [El The lines indicate s pectra for models ex plaini ng the knee due to the max i mum energy attained during the accel- 
erat ion process according tolSveshnikoval ( 2003ft ( — ) . iBerezhko and Ksenofontov ( 1999ft ( — ), Stanev et al. ( 1993ft (• • •), 



and iKobavakawa et al. ( 20021 ) f — ) 



Fig-HU The lines indicate s pectra for mode l s expl ainin g the knee as effect of leaka ge fro m the Galaxy d uring the propa- 
gation process accord i ng to Horandel et al. ( 20071 ) ( — l. lQgio and Kakimotol ( 2003ft ( — ). lRouletJ ( 2004ft (• • •), as well as 



Volk and Zirakashvilil (j2003h (— ■). 



three between the different interpretations. However, it is 
evident that the knee in the all-particle spectrum is caused 
by a depression of the flux of light elements. The KAS- 
CAPE results are illustrated in Figs. fT3l and fT4l 

In the figures also results from further air shower experi- 
ments are shown: EAS-TOP and GRAPES-3 derived spec- 
tra from the simultaneous observation of the electromag- 
netic and muonic components. HEGRA used an imaging 
Cere nkov telescope sy s tem t o measure the primary proton 
flux dAharonian et all Il999l Spectra for protons and he- 
lium nuclei are obtained fro m emulsion chambe rs exposed 
at Mts. Fuji and Kanbala dHuang et all l2003ft . The Ti- 
bet group performs measurements with a burst detector as 
well as with emulsion chamber s and an air shower array 
( Amenomori et alll2000bll2003h . The HESS Cerenkov tele- 
scope system derived for the first time an energy spectrum 



measuring direct Cerenkov light dAharonian et all 120071 ) 



The idea is to register the Cerenkov light of the primary 



nuclei before the air shower starts dKieda et alll200ll ). Re- 
sults for iron nuclei are shown. 

Pespite of the statistical and systematical error of the in- 
dividual experiments, over the wide energy range depicted, 
an overall picture emerges: the measurements seem to fol- 
low power laws over several decades in energy with a break 
at high energies. The bends of the spectra at low energies 
is due to modulation in the hcliosphcric magnetic fields. 

Frequently, the question arises whether the energy spec- 
tra of protons and helium have the same spectral index. 
Puc to spallation of nuclei during their propagation and 
the dependence of the interaction cross section on A 2 / 3 one 
would expect a slightly flatter spec trum for helium nucle i 
as compared to protons. Following iHorandel et al.l (|2007f ) 
the difference should be of the order of A7 « 0.02. How- 
ever, a fit to the experimental data yields 7 P = —2.71 ±0.02 
and 7ff R = -2.6 4 ± 0.02 (grey lines in Figs.Ql and fT4| 
( Horandel . 2003aft . yielding a diference A7 — 0.07. 
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The spectra shown for carbon, oxygen, and the CNO 
group have a threefold structure in resolution. At the low- 
est energies, clean spectra could be resolved for carbon and 
oxygen only (single-element resolution). At intermediate 
energies, the data represent the CNO group. The "CNO" 
group derived from air shower experiments contains prob- 
ably more elements than just CNO due to the limited reso- 
lution in In A. Thus, there is a slight step in the flux at the 
transition from direct to indirect measurements. 

The results shown for iron are single element spectra (as 
obtained e.g. by the ATIC, CRN, HEAO-3, and TRACER 
experiments) as well as iron group data. It may be worth 
to draw special attention to the behavior of the spectra de- 
rived with the model QGSJET from the KASCADE and 
GRAPES-3 data. Both yield a depression of the flux at 
10 6 GeV. Since both experiments obtain a similar effect 
there might be some anomaly in the electron muon corre- 
lation for QGSJET around 1 PeV. It may be noted as well 
that for both experiments QGSJET favors a lighter com- 
position with respect to the SIBYLL results (see Fig.[9|), 
i.e. an interpretation of the measured electron-muon distri- 
butions with QGSJET yields an increased flux of protons 
and helium nuclei, while on the other hand, the CNO and 
iron groups are found to be less abundant in comparison to 
SIBYLL (see Fig.UHJ. 

For oxygen and iron it should be pointed out that the 
spectra of the TRACER experiment (with single-element 
resolution) almost reach the energy range of air shower 
measurements. It is important to realize that with a next- 
generation TRACER experiment overlap in energy between 
direct measurements with single-element resolution and air 
shower arrays could be reached, which is important for an 
improvement of the understanding of hadronic interactions 
in the atmosphere. 

The spectra according to the poly-gonato model are indi- 
cated in the figures as solid grey lines. It can be recognized 
that the measured values are compatible with breaks at en- 
ergies proportional to the nuclear charge Ez — Z ■ 4.5 PeV. 

For the proton component it can be realized that the 
spectrum exhibits a relatively pronounced knee. The spec- 
tral index changes by A7 pa 2.1 from 71 = —2.71 at low en- 
ergies to 72 = —4.8 above the knee. A similarly strong knee 
is visible for the helium component. For the iron component 
no clear cut-off is visible yet. An improvement of the sit- 
uation is expected by the KASCADE-Grande experiment 
( Navarra et al.l . 120041) . which is expected to unravel the en- 
ergy spectra for elemental groups up to about 10 18 eV. 

It is of particular interest to compare the measured spec- 
tra to various astrophysical scenarios. Spectra for four mod- 
els, describing the acceleration of cosmic rays in supernova 
remnants are shown in Fig. 1131 In these scenarios the max- 
imum energy attained during the acceleration process is 
responsible for the knee in the energy spectrum. The mod- 
els make different assumptions on the properties of super- 
nova remnants and, consequently, yield different spectra at 
Earth. 

The calculations by iBerezhko and Ksenofontov (1999) 



are based on the nonlinear kinetic theory of cosmic-ray ac- 
celeration in supernova remnants. The mechanical energy 
released in a supernova explosion is found in the kinetic en- 
ergy of the expanding shell of ejected matter. The cosmic- 
ray acceleration is a very efficient process and more than 
20% of this energy is transfered to ionized particles. The 
resulting spectra at Earth are represented as dashed lines 
in the figure. 

A threefold o r igin o f energetic cosmic rays is proposed 
bv lStanev et al. ( 19931 ). In their mod el particles are ac cel- 
erated at three different main sites ( BiermannL 19931 ): 1) 
The explosions of normal supernovae into an approximately 
homogeneous interstellar medium drive blast waves which 
can accelerate protons to about 10 5 GeV. Particles are ac- 
celerated continously during the expansion of the spheri- 
cal shock-wave, with the highest particle energy reached 
at the end of the Sedov phase. 2) Explosions of stars into 
their former stellar wind, like that of Wolf Rayet stars, ac- 
celerate particles to higher energies. The maximum energy 
attained depends linearly on the magnetic field and maxi- 
mum energies E max = 9 ■ 10 7 GeV for protons and E max = 
3 • 10 9 GeV for iron nuclei are reached. 3) For energies ex- 
ceeding 10 8 GeV an extragalactic component is introduced 
by the authors. The hot spots of Fanaroff Riley class II ra- 
dio galaxies are assumed to produce particles with energies 
up to 10 11 GeV. The spectra are indicated as dotted lines 
in Fig.EH 

A slightly modified version of the diffusive accelera- 
tion of particles in super nova remnants is considered by 



Kobavakawa et al. (|2002h . Standard first order Fermi ac- 
celeration in supernova remnants — with the shock nor- 
mal being perpendicular to the magnetic field lines — is 
extended for magnetic fields with arbitrary angles to the 
velocity of the shock front. The basic idea is that particles 
are accelerated to larger energies in oblique shocks as com- 
pared to parallel shocks. The spectra obtained are plotted 
as dash dotted lines in the figure. 

A usual way to increase E max is t o enlarge the mag- 



netic field B, see e.g. the approach by Kobavakawa et al 



(2OO2]). On the other hand, as consequence of recent ob- 
servations also the parameters of the supern ova explosion 
itself can be varied. This is the basic idea of ISveshnikoval 
(|2003l ). to draw up a scenario in which the maximum en- 
ergy reached in SNR acceleration is in the knee-region of 
the cosmic-ray spectrum (« 4 PeV), using only the stan- 
dard model of cosmic-ray acceleration and the latest data 
on supernovae explosions. Based on recent observations the 
distribution of explosion energies and their rates of occur- 
rence in the Galaxy are estimated. The observed spectrum 
in the Galaxy is obtained as sum over all different types 
of supernova explosions, integrated over the distribution of 
explosion energies within each supernova group. The corre- 
sponding spectra are shown in Fig.[l3]as black solid curves. 
The structures seen in the energy spectra are a consequence 
of different supernova populations. 

Typically, in such models the composition is normalized 
to direct measurements at energies around 10 3 to 10 4 GeV, 



1G 



which explains why the models predict very similar fluxes in 
this region. On the other hand, approaching the knee, quite 
some differences are visible. The least prono unced change 
in the spectrum is obtained in the approach bv lStanev et al 



( 1993f ). while the most dr astic change in the spectra is ob- 
tained in the approach bv lSveshniko va (2003). The predic- 
tions for the mean logarithmic mass of these models are 
summarized in the upper panels of Figs.[9land ll2l 

Examples for a second group of scenarios, in which the 
knee is caused by leakage of particles out of the Galaxy, are 
c ompiled in Fi g. [Ml 



Qgio and Kakimotol ( 2003f ) consider a regular magnetic 



field in the Galaxy following the direction of the spiral arms. 
In addition, irregularities of roughly the same strength are 
supposed to exist. It is assumed that the regular and irreg- 
ular components have about the same field strength B ks 
3 /iG and that both decrease exponentially with a scale 
hight of 1 kpc. The scale length of the irregularities is es- 
timated to about Li rr « 50 pc. The spectra obtained are 
illustrated in Fig. [14] as dashed lines. 

Similar to the model discussed previously, iRoulefJ (|2004 
considers the drift and diffusion of cosmic-ray particles in 
the regular and irregular components of the galactic mag- 
netic field. Again, a three-component structure of the mag- 
netic field is assumed. The regular component is aligned 
with the spiral arms, reversing its directions between con- 
secutive arms. This field (with strength Bq) will cause par- 
ticles with charge Z to describe helical trajectories with a 
Larmor radius Rl — p/(ZeBo). Secondly, a random com- 
ponent is assumed. This will lead to a random walk and dif- 
fusion along the magnetic field direction, characterized by a 
diffusion coefficient Dh oc E m . The diffusion orthogonal to 
the regular magnetic field is typically much slower, however, 
the energy dependence of D± is similar to Du. The third 
component is the antisymmetric or Hall diffusion, which is 
associated with the drift of cosmic rays moving across the 
regular magnetic field. The antisymmetric diffusion coef- 
ficient is Da ~ tlc/3 cx E. For the calculations a source 
spectrum dQ/dE oc E~ aa , with a constant index for all 
species a s = 2.3 is assumed. Below the knee, where trans- 
verse turbulent diffusion dominates, the observed spectral 
index will be a « a s + 1/3, while in the drift dominated 
region above E ■ Ek, a ~ a s + 1 is obtained. This results 
in spectra shown in Fig.[JJ]as dotted curves. The change of 
the slope for individual components Aa ~ 2/3 is relatively 
soft. 

A similar approach is followed bv lHorandel et all ( 20071 ). 
This scenario is, as the two model s discussed above, based 
on an idea bv lPtuskin et all (|l993l ). Solutions of a diffusion 
model are combined with numerically calculated trajecto- 
ries of particles. While in many models the composition is 
normalized to values observed at Earth, in this approach 
the composition at the source is assumed to be equal to 
the abundance of elements in the solar system weighted 
with Z 3 ' 2 . This choice is arbitrary to a certain extent, but 
may be motivated by a higher efficiency in the injection 
or acceleration processes for nuclei with high charge num- 



bers. From the calculations follows a relatively weak depen- 
dence of the escape path length on energy oc E~ ' 2 , even 
weaker a s the one obtained in the model discussed above 
(jRouletl . l2004h . This necessitates relatively steep spectra 
at the sources oc E~ 2 - 5 . However, taking reacceleration of 
particles in the Gala xy into account would lead to flatter 
spectra at the source (Ptuskin, 20061 ). being again compat- 
ible with the standard Fermi picture and the TeV-7 ray ob- 
servations. The spectra obtained are represented in Fig. 1141 
as black solid lines. 

Reacceleration of c osmic-ray particles in the g alactic 
wind is discussed by IVolk and Zirakashvilil ( 20031) . The 
wind is mainly driven by cosmic rays and hot gas generated 
in the disk. It reaches supersonic speeds at about 20 kpc 
above the disk, and is assumed to be very extended (sev- 
eral 100 kpc) before it ends in a termination shock. Due 
to galactic rotation the differences in flow speed will lead 
to strong internal wind compressions, bounded by smooth 
cosmic-ray shocks. These shocks are assumed to reacceler- 
ate the most energetic particles from the disk by about two 
orders of magnitude in rigidity, ensuring a continuation 
of the energy spectrum beyond the knee up to the ankle. 
A fraction of the reaccelerated particles will return to the 
disk, filling a region around the galactic plane (several 
tens of kpc thick) rather uniformly and isotropically. A 
maximum energy E max w Z ■ 10 17 eV is obtained and the 
authors conclude that the knee in the all-particle spectrum 
cannot be the result of the propagation process, instead it 
is supposed to be a feature of the source spectrum itself. 
It is pointed out that it is possible to explain the continu- 
ation of the cosmic-ray spectrum above the knee up to the 
ankle in a natural way, by considering the dynamics of the 
interstellar medium of the Galaxy and its selfconsistent 
extension into a large-scale halo by the galactic wind. The 
authors conclude further that within this picture there is 
no way to produce higher energy cosmic rays, their sources 
must be of a different nature. 

Comparing the spectra of the differ ent models shown in 
Fig. [Til it is seen that the approach bv lQgio and Kakimoto 
(|2003f ) yields the weakest change of the spectra in t he knee 
region . Mo re pronounced change s are obtained bv iRouletl 
(|2004[ ) and lHorandel et alJ (|2007h . The predictions for the 
mean logarithmic mass of these models are summarized in 
the lower panels of Figs.l9"land[T2"l Further models, proposed 
to expl ain the knee and their predictions are discussed else- 
where (|Horandell 12004 120051 ) . 

In general, it may be remarked, that the relatively strong 
change in the measured spectra (at least for protons and 
helium) is not reproduced by the theoretical approaches 
shown in Figs.[T2] and [JJ] It seems to be difficult to re- 
produced such a relatively pronounce d knee as seen in the 
measurements. iHorandel et alJ (|2007l ) tried to model the 
relatively sharp knee by a combination of the maximum en- 
ergy attained in the acceleration process and leakage from 
the Galaxy. This seems to be quite promising, as can be 
inferred from the figure. This, indeed, maybe a hint that 
the knee is caused by a combination of two effects, namely, 
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Fig. 15. Flux of cosmic rays at high energies to illustrate the transi- 
tion of galactic to extragalactic cosmic rays. The data point s are av - 
erage experimental values of the all-particle flux (Horandcl, 2003a). 
The black lines represent possib le fluxes of extra galactic cosmic rays 

according to mod els by iBerezinskv et al.l J2004D ( , two extreme 

cases are shown), Iffillasl ||2005| ) (■ ■ ■) lAtovan and Dermerl d2006l) 
( — ), and the poly-gonato model ( — ) ijHorandel l2003al ) . The grey 
lines give possible contribu tions of galactic cosmic rays according to 
lAtovan and Dermerl l|2006h and the poly-gonato model for elements 
with Z = 1 - 28 and Z = 1 - 92. 

the maximum energy reached during SNR acceleration and 
leakage from the Galaxy. 

6. Transition to Extragalactic Cosmic Rays 

All models for the origin of the knee have in common that 
they predict spectra for individual elements of galactic cos- 
mic rays which exhibit breaks proportional to the charge 
or mass of each element. This ultimately leads to the fact 
that above a certain energy no more particles exist. On the 
other hand, the measured all-particle flux extends up to at 
least 10 20 eV, and the highest-energy particles are usually 
being considered of extragalactic origin. The Larmor ra- 
dius of a proton with an energy of 10 20 eV in the galactic 
magnetic field is w 36 kpc, comparable to the diame- 
ter of the Galaxy. This emphasizes that such high-energy 
particles are of extragalactic origin. The transition region 
from galactic to extragalactic cosmic rays is of particular 
interest, key features are the origin of the second knee and 
the ankle. 

Different scenarios for the transition from galactic to ex- 
tragalactic cosmic rays are discussed in the literature. Some 
recent ideas for a possible flux of the extragalactic compo- 
nent are summarized in Fig. 1151 In these models a possible 
contribution of galactic cosmic rays is estimated and sub- 
tracted from the observed all-particle flux. The different 
assumptions on the behavior of the tail of the galactic flux 
to highest energies yield different estimates for an extra- 
galactic component. 



Reviewing the properties of cosmic rays accelerated in 
supernova remnants, Hillas finds that a conservative esti- 
mate of the maximum energy achieved during the accel- 
eration in supernova remnants is not sufficient to explain 
the all-particle flux up to 10 17 eV. He introduces a sec- 
ond (galactic) component, accelerated by some type II su- 
pernovae, to explain t he observed flux at energies above 
10 16 eV (jHillasl . l2005h . He finds a contribution of extra- 
galactic particles (with mixed composition) as illustrated 
by the dotted line in Fig.[T5l 

In the poly-gonato model a significant contribution of 
ultra-heavy elements (heavier than iron) to the all-particle 
flux is proposed at energies around 400 PeV (jHorandel . 
2003al iHorandel et all [2007h . In this approach the second 



knee is caused by the fall-off of the heaviest elements with 
Z up to 92. It is remarkable that the second knee occurs 
at Ei n d ?s 92 x Ek , the latter being the energy of the first 
knee. The flux of galactic cosmic rays according to the poly- 
gonato model is shown in Fig. [TS] as grey solid lines for ele- 
ments with Z = 1 — 28 (conservative approach) and for all 
elements [Z = 1 — 92). To explain the observed all-particle 
flux, assuming a strong contribution of galactic ultra-heavy 
elements, an extragalactic component, as illustrated by the 
solid line is required. The galactic and extragalactic con- 
tributions according to this scenario are shown as well in 
Fig. [8] as solid (galactic) and dotted (extragalactic) curves. 

The dip seen in the spectrum between 10 18 and 10 19 eV, 
see Fig. [51 is proposed to be caused by electron-positron 
pair production of cosmic rays on cosm ic microwave back - 
ground photons p + 73^ — > p + e + + e~ ([Berezinskvl I2005T ) . 
The flux of galactic cosmic rays is extrapolated from mea- 
surements of the KASC ADE experiment. Two ex treme sce- 
narios are considered by IBerezinskv et"aL (2004). yielding 
extreme cases for a possible extragalactic flux, indicated in 
the figure by two dashed curves. 

An alternative scenario is proposed by Atoyan and Dermer 
( 20061 ). Supposedly, relativistic shocks from a recent (« 
1 Myr old) galactic gamma ray burst (« 1 kpc distance 
to Earth) accelerate galactic cosmic rays up to energies of 
the second knee. At higher energies extragalactic gamma 
ray bursts are proposed to sustain the observed all-particle 
flux. The corresponding galactic and extragalactic compo- 
nents are indicated in the figure as dash dotted grey and 
black lines. 

As seen in Fig.[l5]the estimates of the various approaches 
are quite similar. The extragalactic component reaches its 
maximum (in this E 2 ^ representation) at energies around 
100 to 400 PeV with the flux values differing by about 
a factor of three. The differences are caused by different 
estimates of the tail of the galactic flux to highest energies. 

New measurements of the mass composition in the energy 
region of the second knee will help to distinguish between 
the different models. 

To reach to energies approaching 10 18 eV the KAS- 
CADE experiment has been enlarged. Covering an area of 
0.5 km 2 , 37 detector stations, containing 10 m 2 of plastic 
scintillators each, have been installed to extend the original 



18 



KASCADE set-up (jNavarra et all 120041 ). Regular mea- 
surements with this new array and the original KASCADE 
detectors, forming the KASC ADE-Grande experimen t, are 
performed since summer 2003 ([Chiavassa et al.1 . 120051 ). The 
objective is to recons truct energy spectra fo r groups of ele- 
ments up to 10 18 eV (jHaungs et al.U2006af ). First analyses 
extend the lateral distributions of e lectrons and muons u p 
to 600 m (jGlasstetter et all . 120051 Iv. Buren etahl . 120051) . 
Based on one year of measurements, already energies close 
to 10 18 eV are reached. 

The Ice Cube/Ice Top experiments are installations 
presently under construction at the South Pole. Ice Cube 
is a cubic kilometer scale Cerenkov detector with the main 
objective to measure very high-energy neutrinos with a 
threshold of a few 100 GeV (|Kesteil2004h . It will consist of 
80 strings, each equipped with 60 digital optical modules, 
deployed in the Antarctic ice at depths of 1450 to 2450 m. 
The ice Cerenkov detector will also register high-energy 
muons (E^ > 300 GeV) from air showers in coincidence 
with a surface array. The latter, the Ice Top experiment, 
is located on the snow surface above Ice Cube, covering 
an area of 1 km 2 . It consists of frozen water tanks, serving 
as ice Cerenkov detectors to detect the electromagnetic 
component of air showers. The detectors form an array 
of 80 detectors on a 125 m triangular grid. It is expected 
that the installations will be completed in the austral sum- 
mer 2010/11 and with this set-up air showers up to EeV 
energies will be investigated. 



7. Summary and Outlook 

In the last decade substantial progress has been achieved 
in the understanding of the origin of (galactic) cosmic rays. 
It has become clear that the knee in the energy spectrum 
at about 4 PeV is caused by a break in the energy spec- 
tra of the light elements and the observed knee is relatively 
sharp. The mean mass of cosmic rays is found to increase 
as function of energy in the knee region. The energy spec- 
tra for individual elements seem to follow power laws over 
a wide energy range with breaks at high energies. The posi- 
tion of the breaks are compatible with a rigidity dependent 
scenario, however, a mass dependent behavior can not be 
excluded completely. Despite of this progress, the exact as- 
trophysical interpretation of the measured data is limited 
by the present understanding of hadronic interactions in 
the atmosphere. 

The measurements exhibit qualitative agreement with 
the "standard picture". That means the bulk of galactic 
cosmic rays is accelerated in shocks of supernova remnants. 
The particles propagate in a diffusive process through the 
Galaxy. However, for an exact quantitative description, 
more detailed models are required. Many contemporary 
models are very detailed in one aspect (e.g. acceleration), 
but crude with respect to all other processes (e.g. propaga- 
tion - and vice versa) . One should aim for integrated mod- 
els, taking into account in detail the injection, acceleration, 



and propagation of the particles. Most likely, the knee in 
the energy spectrum is caused by a combination of the max- 
imum energy reached during acceleration and leakage from 
the Galaxy during propagation. However, at present, some 
more exotic models cannot be excluded completely. 

Within the next decade, new measurements in the en- 
ergy region around the second knee should be able to clarify 
the mass composition of cosmic rays in this region. These 
results are expected to contribute to the understanding of 
the end of the galactic cosmic-ray spectrum and the tran- 
sition to extragalactic cosmic rays. 
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